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INTRODUCTION 


It is universally assumed in the text-books on 
suspension bridges that the shape of the curve 
assumed by the cable isa farabola on the suppo- 
sition that the weight of the structure and the loads 
are uniformly distributed in a horizontal direction. 
Neglecting for a moment the moving loads, the 
weight of the structure is made up of the weight of 
the cable and the suspenders, and the weight of the 
roadway and the stiffening strusses, it being assumed 
that the stiffening is not of the braced-chain type. 
We could simply classify this as the weight of the 
cable and the weight of the roadway, the former 
being uniformly distributed along its length and the 
latter horizontally. It is well-known in statics that 
if the former be inappreciable and can be neglected 
the curve of equilibrium is a parabola; while if 
the latter be neglected in comparison with the 
cable weight the curve is a catenary. It is hardly 
necessary to point out that the curve will still 
be a parabola even if the former be not neglected 
but taken to act horizontally. I propose to examine 
in this paper how far this assumption is justifiable 
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and what is the error made by the assumption in 
calculations relating to suspension bridges. In view 
of the fact that every book, either on applied 
mechanics or suspension bridges, makes this 
assumption without any attempt at Justification, it 
is worthwhile to investigate how far it is correct. 
Before proceeding to this discussion, it might be 
observed that greater errors are likely to result from 
this assumption, the greater the weight of the cable. 
The following *table giving a comparative statement 








Net weight of | Gross fixed 
Name. Span in ft. chain per ft. load per ft. 
run. | __ of span. 

Clifton 223 1560 3171 
Freiburg 870 167 760 
Niagara s213 820 2032 











he 


for three suspension bridges shows that the cable 
weight is quite appreciable when compared with the 
gross fixed load and it may consequently be of some 
practical use to undertake this investigation. 
Further, it is also evident that large values of the 
ratio of sag to span introduce greater errors in the 
assumption and in practice this ratio varies from 
0.2 to 0.05 such that it might be as high as 1/5. 


\ 1. THE GENERAL DIFFERENTIAL EQUATION 


To determine the form of the cable for any 
loading, consider vertical loads applied on a cable 
suspended between two points. The cable will 
assume a definite polygonal form determined by the 





“ Table from Rankine : Civil Engineering. p. 582, 
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relations between the loads. The end components will 
be inclined and will have horizontal components H., 
Simple considerations of static equilibrium show 
that H will be the same for both end reactions and 
will also equal the horizontal component of the 
tension in the cable at any point. H is called the 
horizontal tension of the cable. It is easily shown 
that the cable curve is the bending moment diagram 
for the applied loads. If the loads are continuously 
distributed the funicular polygon becomes a conti- 
nuous curve and if w be the load per horizontal 
linear unit at any point having the abscissa x, we 
can obtain the equation for the vertical shear at any 
section of the span and on differentiating this 
equation we obtain, 

d’y 

dx” 





— — w 


and on taking the y-axis vertically upwards this 
becomes 


ay 
dx* 
which might be taken as the differential equation™ of 
the equilibrium curve in the general case. Let g, be 
the weight of the cable per unit length (supposed 
uniform) and w. the total uniform moving and dead 
loads per horizontal foot ; equation (1) then becomes 
d’y 

a dx? 

where « is the inclination of the cable to the horizon- 
tal at any point. Egn. (2) is the general differential 
equation of the problem. If go be assumed to be 


H =w (1) 


= g. secX + Wo (2) 








* See for e.g., W. H. BuRR: Suspension Bridges p. 61. 
Eqn. (13) or JOHNSON, BRYAN and TURNEAURE: Modern 
Framed Structures. "Part II, pp. 187-88. We shall use 
throughout the notation of the former book, 
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uniform horizontally instead of along the curve, we 
can reduce the equation to 


2 
H a = w (w =—Wo a Bo) 


and taking the initial conditions at the lowest point 
of the cable to be: 


dy 


n=O; - =0, y = 0, this equation gives 
X 
nye ave 
WX WX 
iy =e a 
eid AO) ee 


which is the standard equation to the parabolic cable. 
The tension T at any point is given by 


T = H Sec « = HA/, f (2) 
x 





If the points of support are at the elevation y == 7 
and the span = J, H takes the value 


a2? na t= Ba/) AE 
oary and T = 7A +(%) 


If L be the entire length of the cable we have 


l/. dy\2]1 1 
2 2 A. 2 " : Kk 
2 f [ 1+ (2) | dx (20+ sin h2u) * 


Developing the radical] by the binomial] theorem to 
three terms and performing the integration indicated 


L 


IEEE te Ee ame Tew dey 
a 5 es 
The equations for y, H, T, Lare the fundamental 
equations for the parabolic cable. 





* in = Fil and wis defined by sin hy = Dts pe 
: , =+n. This el 
form for L is to be found in Steinman: Suspension Brides ier 
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If in equation (2) we neglect w, we obtain 


2 
HEN=¢, sec & =A/ 14 es eC Y) 
6 


Fi hie 


This is the differential equation of a catenary and 
its solution gives 


1 Cx tot : 
y= ( ny a where c= g./H 
and we can derive corresponding equations for H, T 
and L as in the parabolic cable. 


I now consider the general differential equation 
(2) where neither w, nor g, can be neglected (it 
must be observed that assuming a uniform horizon- 
tal distribution of g, is mathematically equivalent to 
neglecting it) and obtain the exact mathematical 
solution. In the next article, I derive an approximate 
form of the solution suitable for practical use and 
enabling us to compare the formule obtained for a 
strictly parabolic cable. 


We shall write (2) in the form 


H®3 2=£o sec © + Wo (replacing « by © ) 


dx? 
d? = +b, = p=” 
or d 2 aq sec So ( a H H ) 


and assume, as always realised in practice, that b>a. 
Putting p= dy/dx, this further reduces to the first 
order differential equation: 


a2 Cl teitb.> (3) 


Since p = tan ©, dp = sec? 6. dO and (3) can be 
written as 
sec? 6.d0 hl he DO ade 
b + asec O mone cos © (a+b cos ©) rae 





a9) as (bd 





a -=ad 
shins a+ bcos 9 ad 
and integrating we have 
ax = log (sec 8 + tan @)- 
b+ 56+ ,/b?—a? sinO 
b (b+a co a/ | 4) 


72. 42 10g . 
rf ae \— 


the constant of integration being given by x = 0 for 6 = 0. 


a+bcos © 





Again vhs a J + p? + bcan be written as 
S 


_ pdp__ =\dvor tan 6 sec? 6 dO 
ay/1+ p? + b aseccO+t+b 





This can be easily integrated and the result is 


2 afte a _atb ]® 
a? 4 a(sec © 1) + blog [ 4b 


taking y = 0 for 8 = 0. 





Equations (4) and (5) give the parametric equations 
of the curve of the cable in terms of 6, but they are 
in quite an unmanageable form for any practical 
purposes. Before deriving approximate equations 
from these, it is instructive to see how these reduce 
to the parabola and catenary for the cases a=O and 
b=0 respectively. 


Case )>= 0: Equation (4) takes the simple form 
ax = log (sec 8 + tan ©) 
and equation (5) reduces to 
ay = sec 6 —1, since the argument of the 


logarithm in equation (5) can never become infinite. 
Eliminating © between these equations we have 
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+ ( ax —ax ) 
2G \ e te 2 


which is the equation to the catenary 


a 


Case a = 0: In this case equations (4) and (5) reduce to 
identities but by dividing them throughout by @ and 
a2 respectively and applying de L’Hospital’s Rule 
for the evaluation of indeterminate forms, they 
take the form 


Seidl 
x=— ft J 
; an © 


1 f 
and y =>, tan’ © ] 


; bx2 ; 
which give y =~, , the equation to the parabola. 


§ 2. APPROXIMATE SOLUTION OF THE GENERAL 
DIFFERENTIAL EQUATION. 


The differential equation (3) is equivalent to the 
two equations 





dy — —— 


dx a/ — 2 


and 5 = a ( 1—z2) + b ( 1—22)3/ 
KE: 





where g = sin ©. Assuming for y and z two power 
series * with undetermined coefficients, using the 
initial conditions x=0, y=0, z = Oand determining 
these coefficients so as to satisfy the above two 
equations, we deduce for y the expression | 


x2 pen" 
y= Lato) © abt (5) 








+ This method is suggested for approximate solutions of 
differential equations in RUNGE UND KONIG: Numeris- 
ches Rechnen. § 87. pp 321—23 
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Asa first approximation we can retain orlly the second 
term in (5), and the equation then illustrates how 
the assumption of non-uniform horizontal distribu- 
tion of g. affects the shape of the parabolic cable. 


I shall now indicate an alternative method of 
deriving equation (5) which is more illuminating, 
Equation (2), viz. 


d 
KW 4p tb (2) 


gives us the parabolic cable when we write @ in 
place of @ v1 + p? which is equivalent to putting 
p= O z¢e., treating the cable as straight. We can 


therefore obtain a better approximation by putting 





for p iy its value as obtained for a parabolic 


cable. From § 1, the equation 


b x2 da 
y= —>— gives p eatin =bx and hence (2) can be 


approximated to 





d 22 
Pe=sa(1+bx )b +b 
Mees ( b2 x2 b' x' 
hej Sat) a ee 
anal Pinal 2 8 are. 


on expanding by the binomial theorem, since bx 
giving the slope at any point can be assumed, for all 
practical purposes, numerically less than unity, and 





the expansion being therefore valid. Retaining 
only two terms, this gives 

d2 gag 

oF = (ato) + 22 


Integrating once and using a =0 if*¥ =0 
x= 

Wy (3 + b)x + ab? x 

dx : 6 
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Another integration, using y = 0 if * = 0 gives 


(a gn b) x? 4 ab? x! 
2 24 


(5) 


which is the same * as equation (5). For use in 


practical calculations it is preferable to express this 
equation in terms of fand/; but we have only one 
condition expressing that (5) passes through the 
points of support ( + 7//,, f)and this does not enable 
the determination of both aand b in terms of these 
quantities, . Let us, however, write abn, so that n 
is the fraction of the cable weight to the horizontal 
load; equation (5) then reduces to 
BMGT SG, x? nb x 
2 24 
Putting x=1/2 & y=f, we have 
b(nt+1)l? nb? l' 
8 384 
e 2 pire 
bi aE (1 ee (6) 


(64 — == 


ee 48. 





f= 


This equation gives a cubic in 6 from which we can 
find b in terms of , f and which might be supposed 
to be known quantities. We can solve this cubic 
equation approximately. 


o 
Writing —= eqn. (6) reduces to 


fies nds? , a cubic ind 


3 [2 
For small values of , a first approximation is 
zh Synth 
fHd&b= a 





* See W. H. Burr: ibid p. 69 where a similar method is 
adopted in a different connection. 


2 


10 


which is a well-known expression in a parabolic cable. 
To obtain a second approximation let us substitute 
this value ) = fin the right hand side and obtain 


< =f—ng (1 is +f )=s (1-n- EL 


+ SF 5s Pan pe 
0nd = (1 -n—t0) (7) 





Also, ay — (n+1) x eb x (8) 
dx 6 


6 being given from (7) in terms of n, f and J. 
Let us now calculate the length of arc of the 


entire cable as modified by the above formule and 
compare the same with the original expression 


y 4 
Larji+ S22) 


obtained in § 1. Denoting the modified length as 
L’ we have 


' L/. | 1 2\4 : ; 
Is =a fiat; 5 (“’) ) a* a being given from (8) 


/ 2 3 va 
= f0-Po+ { (a5) et i Jie 


Expanding by the binomial theorem and retaining 
terms up to x‘, this reduces to 
dy 


F he + h )2x2 +b) 4qh2~ 8 
ie a2 ff + tet oye { ») tab? — 3(a +b) } x 
=] [: ais (a +b)? | (a+b) { Aa =3 (a +6) } “| 


24 
24 1920 











1] 


on performing the integration indicated. Now, to 
2» . , 
substitute a= bn and b= f ta) in this ex- 


pression, we can write 
2 2\2 
ben Oth es 
l Se 


64 f2 Sif 2 ia 
=" (ian — Sas , retaining only first power of n, 





imi WSR ioe 16 nf 2 
ithiniatiyi 6 = =44— : 


474 
rL'= b2 12 roa 
Now ifs = (1 + 21) oan 


(3n°+3 n —n+)3 | 


Po b2]° mae 
ct {i+ at 2n) oon (3+2n)} 


neglecting terms in 72, 1°, etc. 





{,, 0212 _b‘14 ne R (1-22) | 
( gf) 1920" 12 so J | 


Substituting in this the values of 62 and 06‘ written 
above, we get finally 


me Bee 32 fieegias fs, 512 ee 
ert a Pai — aah alan, aaa 





Hence L’ - L=/ = n f4/14 (neglecting the last term) 
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In practice it is found that all ratios of f/t will 
be between .2 and .05 and the following table * 
gives values of L/l corresponding to different values 


of f/1. 














fil LH 
af 1.09643 
ly 1.07566 
sL9 | 1.05676 
AZ5 1.040104 
2 | 1.026027 
075 1.014797 
05 | 1.006627 
* Vide Burr : ibid. (p. 73) 
‘The above expression is = n. Z; 


gives us an indication of the error committed in 
treating go as a uniform horizontal load. Taking from 
the table in the Introduction the case of the Freiburg 
suspension bridge for which »=+ (approximately) let 
us calculate the deviation for the case #//=.2; we 
easily find 
= ="005262 

such that the corrected value would be 1.10169, 
showing that although the modification required is 
not great, it is of the same order as the difference 
between the values of L/7 calculated as in the above 
table and the value calculated from the formula + 
exact for a parabolic cable, viz. 


1 : 
Lil= > (20 + Sin h2u ) 








t Vide Burr: ibid. (p. 7) 


1.3 
§ 3. OTHER DIFFERENTIAL EQUATIONS 


In the previous investigation we have considered 
the weight of the cable alone and neglected that of 
the suspenders. Let g’ be the weight of the suspenders 
per square unit of vertical plane of the cable so that 
g’y is weight per linear horizontal foot of that part of 
the suspenders above the vertex of the cable. The 
part of the suspenders below the vertex forms a 
uniform load over the span and might beassumed to 
have been included in W,. The general differential 
equation becomes 





d’y 
ea it pitotey 





: d 
Ai oo nee V1 +p), + .b icy (c=g'/H ) 


Making # and y the dependent and independent 
variables respectively, this can be reduced to the first 
order differential equation. 





eet V1 +p°+ b+ cy 


of which an exact solution is complicated. It can 

however be solved by approximate methods by 

substituting for y and # the expressions as given 
dx 


by the parabolic cable. 


2 
Putting y="* and Sor, equation (10) reduces to 





ay — Nips a+ bb gh 
di? a a/ 1+ b2 x D 

b2x2 bex 2 

= 2 +b: 
a ( Bs 2 ) 2 


=(a+s)+ x(b +) 
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Successive integrations giv@& — 





dy (ath) go Oe 
dx 6 
2 4 
a (a+b)x 4 ,(abte) x 
Zz 24 


which expresses the deviation from the parabolic 
form. 


I now consider, finally the case where the cable 
is not necessarily uniform. Let g, be the weight per 
ft. run at the vertex and g, at any point where 
inclination of the cable to the horizontal is 6, 
supposing that the normal section of the cable at 
any point is proportional to its total stress, we have 


& x=8 0 SecO xX SecO =g,_ Sec? O. 
and the differential equation becomes 


d2 
“* =a SectO +b+cy (11) 
dx? 

Putting 62 ="¥and y= <f x? 


this can be written 
rite 
d Y=(a+b)+ (AF 4he% 2 


dif i ie 
Ss oe ; 
giving y= (a+b) +kx 
2 
where k= -L (ee tof 
l 12 
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An exact solution of equation (11) is a matter of 
some difficulty. 


I hope to return in a further paper to more 
detailed discussions of equations (10) and (11) and 
also discuss, if any, the bearing of the results of § 2 
on the problems of stiffening of the structure and the 
deflection of the bridge system, 


